uem*rr 

TECHNICAL REPORT S F C T fQM 

NA¥AL ^o^graduai l school 
MC^WT'E^ tTY. CALIFOHNfA 93940 



NPS-61Sd76071 



NAVAL 



POSTGRADUATE SCHOOL 

Monterey, California 




AN INTRODUCTION TO THE SONAR EQUATIONS 
WITH APPLICATIONS 

A. B. Coppens and J. V. Sanders 

July 1976 
NPS Publication 

Approved for public release; distribution unlimited 

FEDDOCS ‘pared for: 

D 208.14/2:NPS-61SD76071 ra.1 Postgraduate School 

iterey, California 93940 




NAVAL POSTGRADUATE SCHOOL 
Monterey, California 



Rear Admiral Isham Linder Jac k R. Borsting 

Superintendent Provost 



The materials presented herein were produced with the support of the 
Department of Physics and Chemistry and the Office of Continuing Education, 
Naval Postgraduate School. 

Reproduction of all or part of this report is authorized. 

This report was prepared by : 



UNCLASSIFIED 

SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 



REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


1. REPORT NUMBER 

NPS-61Sd76071 


2. GOVT ACCESSION NO. 


3. RECIPIENT'S catalog number 


4. TITLE (and Subtitle) 

AN INTRODUCTION TO THE SONAR EQUATIONS WITH 
APPLICATIONS 


5. TYPE OF REPORT 6 PERIOD COVERED 

NPS Publication 


6. PERFORMING ORG. REPORT NUMBER 


7. AUTHOR(s) 

Alan B. Coppens and James V. Sanders 


8. CONTRACT OR GRANT NUMBERS 


9. PERFORMING ORGANIZATION NAME AND ADDRESS 

Naval Postgraduate School 

Monterey, California 93940 Code 61Sd/Cz 


10. PROGRAM ELEMENT, PROJECT, TASK 
AREA 6 WORK UNIT NUMBERS 


1 1. CONTROLLING OFFICE NAME AND ADDRESS 

Naval Postgraduate School, Code 61 
Monterey, California 93940 


>2. REPORT DATE 

July 1976 


13. NUMBER OF PAGES 
121 


14. MONITORING AGENCY NAME a ADDRESS^// different from Controlling Office) 


IS. SECURITY CLASS, (of this report) 

UNCLASSIFIED 


15«. DECLASSI FI CATION/ DOWN GRADING 
SCHEDULE 


16. DISTRIBUTION STATEMENT (of thia Report) 

Approved for public release; distribution unlimited. 


17. DISTRIBUTION STATEMENT (of the abetract entered in Block 20, if different from Report) 


10. SUPPLEMENTARY NOTES 


19. KEY WORDS (Continue on reverse aide it necessary and identify by block number) 

Sonar Attenuation 

Underwater Acoustics Detection 

Sonar Equations 
Underwater Sound Propagation 
Transmission Loss 


20. ABSTRACT (Continue on reverse side it necessary and identify by block number) 

This report provides an introduction to the SONAR equations for those 
interested in underwater sound as applied to ASW but lacking either the 
mathematical background or the time for a more rigorous presentation. Earliei 
versions of these notes were developed for Continuing Education courses 
presented at Moffett Field, California, and Naval Torpedo Station, Washington. 
Additionally, these notes have been in demand for certain courses at the 
Naval Postgraduate School. While this is the text for these courses and 



DD 1 JAN*73 1473 EDITION OF t NOV 65 IS OBSOLETE 

S/N 0102-014- 6601 | 



UNCLASSIFIED 



SECURITY CLASSIFICATION OF THIS PAGE Dmta Entered) 



UNCLASSIFIED 

^LCURITY CLASSIFICATION OF THIS PAGEW>®n Dmtm Entmrmd) 



20. should be supplemented by lectures, we have attempted to design the 
material so that it is reasonably self-explanatory, communicating many of 
the essential concepts without requiring extensive verbal amplification. The 
unusual format has been deliberately chosen to facilitate these goals, and 
our experiences in presenting these materials have seemed to justify this 
choice. It is assumed that the reader has some familiarity with trigonometric 
functions and either has or will develop with the aid of the appendix the 
facility of handling scientific notation and logarithmic operations. 




UNCLASSFIED 



SECURITY CLASSIFICATION OF THIS PAGEfWhen Dmtm Entmrmd) 



TABLE OF CONTENTS 



I. Elements of Acoustical Propagation in the Ocean 
Fundamental Acoustical Quantities 
Sound Pressure Level 
Transmission Loss 
Speed of Sound in the Ocean 
Ray Acoustics 

II. Transmission Loss Models 
The Mixed Layer 
Surface Interference 
The Convergence Zone 
Bottom Bounce 

The Reliable Acoustical Path 

The SOFAR Channel 

Acoustical Normal Modes 

Rays, Normal Modes, and Reality 

Summary 

III. Combination of Acoustical Signals 

IV. The SONAR Equation 
List of Symbols 
Passive Sonar 

Source, Noise, and Masking Levels 
Array Gain and Directivity Index 
The Passive SONAR Equations and FOM 
Detection Threshold 
Doppler Shift 
Examples 

Active SONAR 

List of Additional Symbols 
Source Level 

Target Strength and Time Stretching 

Masking Level (Noise and Reverberation) 

Doppler Shift 

Detection Threshold 

Summary 

Examples 

V. Appendix A: Logs 

Appendix B: Index of Useful Tables and Formulas 



Page 

3 

11 

13 

21 

25 

29 

34 

40 

45 

47 

48 
50 

54 

55 

56 

65 

66 
67 



93 



110 

121 



1 



PURPOSE 



One example of a SONAR equation, for the passive detection of a target 
which is radiating sound into the ocean, is 

I 

TL >, NSL + 10 tog w - DI + DT 

Source Level of the target being detected passively 
Transmission Loss as the signal propagates to the detector 
Noise Spectrum Level of the ambient noise in the ocean 
Frequency Bandwidth of the detecting system 
Directivity Index of the detecting system 
Detection Threshold of the detecting system for a specified 
probability of detection and a specified false-alarm rate 

The purpose of this course is to provide the physical and conceptual 
background necessary to understand the meaning of each of these terms 
and how they are combined to form the SONAR equations. At the 
end of this course, the student should be able to apply the appropriate 
SONAR equation to a given problem and estimate such things as maximum 
detection range, optimum receiver depth, necessary bandwidth, etc. 



EL - 
where 
EL = 
TL = 
NSL = 
w = 
DI = 
DT = 
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"Totn I P* e s s v r e (Pa) 



ACOUSTIC PRESSURE 







Acp^sifc t - (P(0 - (P *>K*) • 

f(-t) tV ie s CmseJ by 
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HYDROPHONE OUTPUT VOLTAGE 





- 'ir(-t) 

Hydrophone sensitivity / x = ( i. ) 

where 'tr(t) j. s the instantaneous voltage at output 
terminal of hydrophone 




h 



hydrophone 



DISPLAYS 




OSCILLOSCOPE 

instantaneous voltage v(t) 




VOLTMETER 

some average voltage 

or 

\ 




BEARING-RANGE RECORDER (active sonar) 

brightness related to the strength of 
the echo returned from each range and 
bearing 



LOFAR-GRAM (passive sonar) 

paper is electrically burnt, with 
darkness related to the strength of 
the signal received at each frequency 
and time 




MONOFREQUENCY SIGNALS 



For ease of discussion, we will at first restrict ourselves to consideration 
of tonal s (monofrequency signals) . Later on, we will see that more compli- 
cated signals can be broken down by FFT, Fourier Analysis, or other filtering 
techniques into a collection of monofrequency signals. 

The time history of the acoustic pressure for a monofrequency signal can be 
represented as 

jc V 

A 

p 

. ^ 

-A 



P = The effective amplitude of the acoustic pressure. This is not the 
same as the peak amplitude A, but is more convenient. We will never 
use the peak amplitude, but it turns out that P = 0.707A for the above 
signal. 

T = The period of the monofrequency signal. Tne period is the time 

interval required for the acoustical signal to go through one complete 
cycle, returning to the configuration it had at the beginning of the 
interval . 




T 




Tf Alt 
>- 
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PLANE WAVES 



jo ~ / A 3//V ^7T^( t £) 

If x = constant (For simplicity, take x = 0.) 

£ - 

A 

A 

A 




If £ = constant (For simplicity, take t = 0.) 







k 



As time increases, the wave moves towards larger x. 



7 



f* 1 h«%y*D 



A WAVS IV SPACE AND TIM 
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Of 



INTENSITY 



"7 “7 




Intensity = power per unit area 

where, "area" is the area of a "window" perpendicular 
to the direction of travel of the wave. 

In terms of measureable quantities, 




where P C 

J 



PC 



= specific acoustic impedance of the fluid 



1.5 x 10 Rayl 
415 Rayl 



for water 
for air 



1 Rayl = 1 kg/m sec 
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SIZES OF THINGS IN WATER 



Speed of sound, frequency, and wavelength: 



The speed of sound in sea water is approximated by 

3 

c ~ 1. 5 x 10 m/sec _+ a few per cent. 

It must be remembered that because of depth, temperature, salinity, and 
contamination there may be variations up to a few per cent in the exact 
value of c. 



Approximate sizes of wavelengths of monofrequency sounds of different 



frequencies are: 



Frequency 
10 Hz = 10 



cps 



50 

100 

500 



1 kHz = 1000 cps 
5 

10 

50 

1 m = 1.09 yd 



Wavelength 
150 m 
30 
15 
3 

1.5 

0.3 

0.15 

0.03 = 3 cm 



Pressure amplitudes: 




2 

1 atmosphere =14.7 lb/ in at the ocean surface 
0 

10 iib (yb = microbar) 



(P 



o 



= lO^yPa (yPa = microPascal) 
increases by about 1 atm for every 33 ft of depth in water. 



Representative acoustic pressures generated near a sonar transmitter: 

P ~ 10^ to 10^ yb or 0.01 to 1 atm. 

If P ~ (f ( , then cavitation is likely and performance nay be degraded. 
Representative weak acoustic pressures detectable by a sonar receiver: 

P 10 - ^ yb or 10 - ^ atm. 
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SOUND PRESSURE LEVEL 



SPL = 20 log(P/P ref ) re P^ f 



0.0002 yb 


(old noise measurements) 


p ref - lub 


(old conventional standard) 


lyPa 


(new Navy standard) 



If P = 1 atm = 10 
194 




If P = 10" 4 atm = 

r 

114 

SPL =<40 
J.40 

If P = 10 -8 atm = 

( 34 
-40 
60 



Examples 

8 yb, then 
dB re 0.0002yb 
dB re lyb 
dB re lyPa 

lO^yb, then 
dB re 0.0002yb 
dB re lyb 
dB re lyPa 

10” ^yb, then 
dB re 0.0002yb 
dB re lyb 
dB re lyPa 
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CONVERSION BETWEEN SPL’s WITH DIFFERENT REFERENCE PRESSURES 



For a given pressure amplitude 



SPL re lyPa = SPL re lyb + 100 dB 
SPL re 0.0002yb = SPL re lyb + 74 dB 



Notice 



0 dB does not mean the absence of sound , but means P = P 

ref 



If 

then 



P = P 



ref ' 



SPL = 20 log <P/P ref ) = 20 log(P ref /P ref ) = 20 log 1 = 0 dB re P^ f 
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TRANSMISSION LOSS 




P(rO 



SPL( r-^ ) 


= 20 log p — — 

ref 




SPL(r 2 ) 


P(r 0 ) 

= 20 log p-i- 

ref 




SPL(r^) 


- SPL(r 2 ) = 20 log 


PCr-jO 


P(r 2 ) 


TL = SPL 


(1 meter) - SPL (r meters) = 



P(l) 

PTrT 
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SAMPLE CALCULATION OF TRANSMISSION LOSS 



If SPL at 1 meter = 


120 dB re lyb 


and SPL at range r 


= 40 dB re lyb 



then TL = 120 - 40 dB = 80 dB 

TL does not depend on the reference pressure. Let’s convert the above SPL's 
to reference lyPa and recalculate the TL 



New SPL at 1 meter = 


220 dB re lyPa 


and SPL at range r = 


140 dB re lyPa 


so that TL = 220 - 


140 = 80 dB 
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TRANSMISSION LOSS FROM SPREADING 



The total transmission loss TL can be considered to arise from two different 
causes : 

1) As the rays of sound propagate out from the source and travel through 
the water, they will bunch together or spread out depending on the properties 
of the speed of sound profile for the water. The transmission loss which arises 
from this effect is termed 

TL = transmission loss from spreading 
Some simple examples of spreading for rays traveling straight lines are 

Spherical spreading 

TL = 20 log r 
g 6 






Cylindrical spreading 
TL = 10 log r 

s 




No spreading 

TL = 0 dB 
g 



In the real ocean, the rays never travel in straight lines, so that these 
above simple examples are unrealistic. We will see, however, that there are 
many situations in which the sound will spread out spherically or cylindrically 
over large portions of its path, and the above simple equations will appear 
in the complete transmission loss equation. 
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SPHERICAL AND CYLINDRICAL WAVES 



SPHERICAL 







CYLINDRICAL 




jb - siv 3 7,J r ( t ~ 

IF ^ » A 



16 



TRANSMISSION LOSS FROM DISSIPATION 
2) There are a variety of other physical effects termed dissipative 
effects, which reduce the SPL of the sounds in addition to the spreading 
of the rays. Since the spreading of rays is described by TL , we will call 
the portion of the transmission loss which results from these other effects 
TL = Transmission Loss from dissipation 

X/ 

There are a number of these which are important in underwater sound propagation: 

(a) Absorption. As sound travels through the water, there is a 
conversion of the acoustical energy into thermal energy because 
of effects which are very similar to friction in moving machinery 
and hysteresis losses in transformers. 

(b) Transmission into the bottom: When a sound beam strikes the ocean 
bottom, part of the beam is reflected back upward in the water 
and another part is sent down into the bottom. The beam which 

is reflected back upward in the water contains less energy 
than the beam which traveled down to the bottom. 

(c) Scattering. The presence of rough boundaries at the ocean surface 
and bottom can lead to a scattering of the sound out of the reflected 
beam. The presence of bubbles and other impurities in the water 
itself can also scatter sound energy out of the path of the ray. 

(d) Leakage. In certain kinds of sound propagation in the ocean, 
like the surface layer, the rays of sound can "leak" out the 
bottom of the layer each time they bend down to graze it. This 
"leakage" results in a loss of energy from the layer. 
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TRANSMISSION LOSS 



In general the total transmission loss TL can be composed of two terms. 



TL = TL 


+ TL. 


g 


1 



where 



TL = transmission loss from geometrical spreading 



TL^ = transmission loss from dissipation 
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SPEED OF SOUND, TEMPERATURE , AND DEPTH 
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These curves show haw the speed of sound is related to depth and 
temperature. The effect of salinity is not included, so these 
curves are appropriate only for water of constant salinity, 
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SPEED OF SOUND FORMULA 






H. (?u T - 0.55 T* 

-f 

I . 3 ( 5 - 3 ?) 

4 

0,017 ^ 



Temperature Correction 
T in °C 



Salinity Correction 
S in parts per thousand 



Depth Correction 
z in meters 



Very accurate formulas for the speed of sound as a function 
of temperature , salinity, and depth are very involved and 
difficult to calculate . The above expression, although 
not highly accurate, is sufficient for many calculations 
and is relatively easy to use* It is accurate to within 
about 0.1 meter per second. 

For a discussion of many of the more precise formulas and 
the discrepancies between them, the interested reader is 
referred to V. A. Del Grosso, J. Acoust. Soc. Am. 56 , 
1084-1091 ( 1974 ). 
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TYPXAl SOUND SPEED PROFILE 
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PROFILE VARIATIONS 



HOURLY 



WEEKLY 



SEASONAL 
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Depth 



Speed of Sound c(z) 



HAYS 



Ran£e r 








RAYS IN CONSTANT GRADIENTS 



e <J> 




/ 

/ 

/ 
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/ 

/ 







^5^ c . + frfr 



t 



is the gradient 



The gradient g can be either positive or negative. 

For a speed of sound profile composed of layers 'within each of 
which the gradient has constant value, then the segment of the 
ray within each layer is the arc of a circle whose center lies 
at the depth z - -c Q /g below the top of the layer (where the 
speed of sound is c Q at the top of the layer and g is the value 
of the gradient for that layer). The radius of the circle is 

C 

- ■ ■ ■■ » 

Cent & 

The ray always bends toward the lower speed of sound. 
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CONSTRUCTION OF RAY PATHS 



Line of centers for bottom layer 

\ 

I ; 

1 

I 

j 

< 

i 
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ACOUSTICAL RECIPROCITY 




If the direction the sound travels over the propagation paths is 
reversed by exchanging the positions of the source and receiver, and the 
propagation conditions remain the same, then the transmission loss 
between source and receiver will remain the same. 




U 

Active 




As a consequence, the transmission loss suffered by the signal in 
reaching the target from an active source will be identical with the 
transmission loss experienced by a signal of the same frequency (and 
duration) radiated from the target when it reaches a receiver. 
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RAYS AND A SURFACE DUCT 



LOSS FROM THE SURFACE DUCT 



When a ray reflects from the ocean surface, some of the energy 
is scattered into different directions and can propagate out of the 
surface layer. This energy is thus lost from the channel, which 
increases the transmission loss. Additionally, energy can diffuse 
out of the bottom of the layer. These effects can be described in 
terms of the "loss per bounce". The number of bounces is the range 
divided by the skip distance r^. 

An empirical formula, resulting from the analysis of a number 
of experiments, has been obtained by W. F. Baker, J. Acoust. Soc. 
am. 57, 1198-1200 (1975)-- 



This equation is valid within a restricted range of parameters, 
given approximately by 



b 




b = loss per bounce in dB 



SS = Sea State 



2 < SS < 5 



3 kHz < f < 8 kHz 



80 ft < layer depth K 200 ft 



30 ft < source depth < 60 ft 
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TRAMS MISSION LOSS IN THE MIXED LAYER 




The depth of the mixed layer is designated by D. 



The depth of the source in the layer is designated by z . 



The transition range r is the distance the sound beam must travel 
until it effectively t ' fills the channel. It is calculated from 



A - (. 05*10 

t 



x / P* 

V D-fc 



(all distances are in meters) 



The skip distance r g is the distance required for the rays which just 
graze the bottom of the layer to complete one full cycle from 

the layer bottom to the surface and back to the bottom. It is 
calculated from 



/■t- = 8.4 »/0*/F 

A ~ 



(all distances are in meters) 



If the loss per bounce is designated by b and the energy loss resulting 
from absorption is designated by a, then the transmission loss for the 
mixed layer can be estimated by tTie following formulas: 



TL * 20 log r + ar within transition range 

TL * 10 log r^ + 10 log r + ar + br/r beyond transition range 



For ranges less than r the sound is spreading out spherically to fill 
the sound channel. For ranges greater than r , however, the sound 

has spread far enough vertically to be reflected ^ from the surface and 
bent back upward from the bottom of the layer. Once this hap ens, then 
the sound is trapped and will spread cylindrically. 
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SAHPLS CALCULATION OF A 



TRALSKtESICl LOSS 



Assumptions: 



Depth of the mixed layer = i *0 n 

Depth of the source * 20 n 

Frecuency cf ti e sound » 2 kHz 

Sea state * 3 

\ e then determine 

Skip distance = 5300 m 

Transition ranpe = 9 u 0 m 

Attenuation constant = 0.00015 dB/rn 

Loss per bounce ■= 3.5 dB/bounce 

The transmission loss formulas are thus 

(a) For ranges less than the transition ran(.e, 

TL = 20 lot r + (1.5 x 10"^)r 

(b) For rar£es greater than the transition rar.pe, 

TL = 10 lop r + ( 8.1 x 10 “k)r + 30 . 
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TRANSMISSION LOSS CURVE FOR SOURCE AND RECEIVER IN TIE DUCT 



Parameters s 

Depth of the mixed layer ■ 1*0 m 

Depth of the source - 20 m 

Frequency of the sound « ? kHz 

Rea state * 3 

Receiver is in the duct 



AW/Fif 

JOO >006 16,000 




. tr 



BUT: These results can be radically changed by effects such as 
Surface Interference (which becomes important at lower 
frequencies and lower sea states), steep gradients below 
the layer, and source and/or receiver near the top or 
bottom of the duct. Additionally, if the frequency of 
the sound is too low, it will not be trapped in the duct. 

As a very rough estimate, the frequency must be greater than about 
(2 x 1C)5 )/d 3/2, where D is in meters and the frequency is in Hz, 
for the sound to be well trapped in the layer and for the above 
model to apply. 
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surface interference 
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Receiver 
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d ~ RECEiysP otrr‘( 

f r /9 /I A//?£- f UO /? / Z 0H7,\L ) 



34 



Acoustic pressure at the receiver 




pressure amplitude ■ A 



A d > A r 



A d * A 

A d » A 



r 

r 



because (l) reflected rays travels a longer path 
(2) scattering from rough surface 
if range large and surface smooth 
if surface very rough 



phase delay ■ fa (due to difference in path length and to phase change 

upon reflection) 
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Transmission loss 



j- f, N smooth surface A d A r 




ran£e (km) 





ranpe (lan) 



Heal sound speed profiles will affect TL only at lar^e 
distances, (hotice the shadow zone predicted at d km.) 



36 



The rar.ge at vhich the received signal is a minimum is given by 



r 



2. h d r 

0-1 C r 



vhere 0 * L ' J , 2 , 3 , ' 



It is possible to use information from the received signal to determine the 
depth, speed, and distance of closest approach of the source. For a source 
moving at constant speed v along, a straight line at a fixed depth h, the 




vhere at t = 0, r * s * the distance of closest approach. 

Combining these equations gives the time at vhich the signal is a minimum as a 
function of frequency. 



i l 1 ( 2 h ol j _£ 

+ W * (-zrr) ? 
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The display that results is shown below 
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Our equation can be rewritten in standard form 



where 




a r 



c S 
PhJ 



I 



and b> ~ 'S/'lr 



This is the equation of a hyperbola 



where a is the value of f at t. 




0 and t/a is the slope of the asymptote. 



If v is calculated from the measured Doppler shifts of the discrete lines, * 
the depth of the source is found from 




and the distance of closest approach from 




Of course, for day-to-day use, nomograms can be developed to simplify the 
calculations. 



* It will be shown later that 

A* ’ 2 f f 

where is the total frequency change of a line. 
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PROPAGATION PATHS FOR S HALL 01' SOURCE Af.D RECEIVER 





(a) Surfade Duct . The isothermal water at the ocean surface forms 
the mixed layer in which sound can be trapped. The continual 
reflection at the rough surface, refraction at the bottom of 
the layer, and leakage provide relatively large energy losses 
so that this path is useful for relatively short ranges. 

(b) Convergence Zone . Vhen the water is deep enough, the sound rays 
which travel down past the SOFAR channel axis can be bent back 
upward before they reach the bottom. These rays, when they 
approach the surface again, tend to be brought together so that 
the sound energy becomes more concentrated, and the transmission 
loss is reduced by a significent amount. Under advantageous 
conditions, these rays will be refocused in second, third, and 
higher convergence zones. 

(c) Bottom Bounce . The reflection of sound off the bottom and back 
to the surface provides a propagation path filling in the ranges 
between the greatest useful range of the mixed layer and the 
innermost edge of the convergence zone. Unfortunately, the 
unpredictable and usually large reflection loss that occurs at 
the bottom causes the transmission loss along this path to be 
rather large. 
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'ransmisrion Loss (dB) 



f - 2 kHz 



a * 0.08 dB/km 




/ /o / oo 1,000 

Range (km) 



(a) Surface duct 

(b) Convergence zone 

( c ) Bottom bounce 



Ref. Urick, Principles of Underwater Sound for Engineers 
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CONVERGE! CE ZONE PARAMETERS 



\<-AC -» 





c Q * the maximum speed of sound above the SOFAR channel axis 
Ac* the difference between c Q and the speed of sound at the elbow 
Z ■ the vertical distance between depths for which c(z) = c Q 
r ■ the ran£e to the conver£ence zone 

\jLk 

Snell's Law reveals that a ray which is horizontal in the upper ocean where 
c(z) ■ c will also be horizontal at a distance Z deeper, and a ray which 
is directed downward where c(z) * c Q must penetrate below the distance Z. 

To insure that enoufh rays are returned to the surface to focus at the 
convergence zone, the total depth of water must exceed Z or Z + D by an 
amount called the depth excess ., d 

If there is no surface layer, then c is associated with the surface 
and Z is measured from the surface. If there is a surface layer, then c q 
is associated with the bottom of the layer and Z measured from the 
bottom of the layer. 



r<— AC 
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C0KVERCS&C5 zone ranch atd transmission loss 



A rough estimate of the range at which a convergence zone forms is given 

fcy 

i 

r C 2 ~ 2.8 Z + Layer effect 



Layer effect ■ 0 if there is no surface duct 

■ r if a surface duct is present and toth source 
aftd receiver are near the top of the layer. 



The accuracy of the above expression depends on how well tue constant- 
gradient approximations fit the true c(z) profile. As might be expected, 
when the elbox-J of the curve is very sharp, and the gradients above and 
below the elbow are relatively constant, then the prediction >111 be 
rather gcod. The width of the convergence zone, over which the rays 
are appreciably focused, is normally about 10 per cent of the range to 
the zone. 

Since the rays begin to spread out spherically as they leave the 
source and are attenuated only by the absoxptive losses in the ocean, 
the transmission loss for convergence zone paths can be written in the 
form 



TL = 20 log r + ar - G 



where C is the convergence zone gain . Typical values of 0 range between 
about 5 and 15 dli when a convergence zone is found. 
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Typical Convergence Zone Langes 



North Pacific 



Surface Temperature-^ 
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33 






* Assumes surface duct absent 
+ Allows for a 200 ftm depth excess 

Ref. Convergence Zone Range Slide Rule, Naval Underwater Systems °enter 



44 



BOTTOM BOUNCE 



For angles of depression beyond about 15°, the radius of curvature of 
a ray is so large that it is a good approximation to assume straight-line 
propagation of the rays over the distances of interest. 




A ray bouncing off the bottom ■will return to the surface at a range 

r •* 2 H tan 0 

Since the distance travelled by the ray is r/(cos 0 ), and since the 
rays spread out spherically and are attenuated only by the absorptive 
losses of the ocean and the reflection at the bottom, an approximate 
transmission loss can be written as 



TL 



20 log 



r 

cos 0 



+ 



ar 

cos 0 



+ BL 



where BL is the bottom loss in dB/bounce. 



I 
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BOTTOM LOSS 



Depan Is 



r/2H 



6 




Depends on bottom composition 

Example: At 2k kHz and 10° grazing anfle 



Mud 

Mud -sand 
Sand-mud 
Sand 
Stony 



16 dB/bounce 
10 
6 
u 
k 



Depends on surface roughness 

Bottom loss goes up as ratio of rou£hness-to-wavelength increases 



Ref. Urick, Princioles of Underwater Sound for Engineers 
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Depth 



RELIABLE PROPAGATION PATH FOR DEEP SOURCE OR RECEIVER 





Reliable Acoustical Path (RAP). 

If the target is shallow and the sensor deep, or vice versa, there 
exists the possibility of a guarenteed acoustical path between the two. 
This path is useful to ranges beyond those attainable by transmission in 
the surface duct. As the deeper of the two (sensor or target) is lowered 
further and further below the layer, the shadow zone indicated in the 
diagram moves out to larger ranges. If there is no mixed layer, then the 
leading edge of the shadow zone is formed by the downward-bending ray 
which just grazes the surface. 

Because the angles for RAP are usually fairly large, the rays can 
be considered to travel approximately in straight lines between the 
target and sensor. For this kind of propagation, the rays spread out 
spherically and the only losses are those from the absorption of sound 
in sea water. An approximate transmission loss formula is then 



TL 



20 log 



r 

cos 9 



♦ ^ , 
cos 0 



where 0 is the angle of elevation of the rays between target and sensor. 

If the negative gradient (below the layer, if one exists) has 
magnitude g, then the range to the shadow zone at the surface is estimated 
by 

i.U V c 0 ( z “ D )/&' ♦ h r s 

where (z - D) is the depth below the bottom of the layer. If there is no 
layer, then D * 0 and there is no skip distance, so r g • 0. 



47 



Depth 



THE SOFAR CHANNEL 




c(z) 



Range 



T 



V 



Just as for the surface duct, the sound spreads out more or less 
spherically until the rays are bent back into the channel, after •which 
the sound spreads cylindrically. Thus, beyond some transition range r^’ 
we should have a spreading term which goes as 10 log r. Those rays 
whose angles of inclination with respect to the SOFAR channel axis are 
large enough that they reflect from the ocean surface or bottom will 
suffer reflection losses and scattering losses and eventually be so 
weak compared to the trapped rays that they can be neglected. Thus, 
for propagation in the channel, there are no losses except absorption, 
so that the transmission loss has the simple form 



The transition range depends on how deep the sound source is in the channel 
and the vertical extent of the channel. For sources placed reasonably 
down in the channel, the transition range is less than the convergence 
zone range r . Notice that the rays cycle up and down in the channel, 
accomplishing^one cycle over an interval very similar to the convergence 
zone range. 



TL = 10 log r + 10 log r^' 



+ ar 



for ranges beyond r^ 1 . 



n 
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TIME STRETCHING IK THE SCFATv CHAM, EL 



time 



L ► 







An interesting fact often observed in the SOFAE channel is that transient 
signals are "stretched" in time as they propagate to large ranges. This arises 
because of the variation in the speed of sound over the vertical extent of the 
channel. For the SCFAR channels commonly encountered in the northern hemisphere, 
the speed of sound increases sufficiently rapidly with distance above or below 
the channel axis that sound traveling oYer those rays which swing over large 
vertical distances, like a and b in the above sketch, reaches the receiver 
at times earlier than sound traveling over the more direct paths c, d, and e. 
Thus, the sound that travels straight down the axis of the channel arrives at 
the receiver last. This leads to the characteristic blossoming of the received 
signal from an explosive charge, the sound first being very weak, and then 
growing in volume until there is a sudden cutoff of sound as the last signals 
come over the slowest paths. 

The exact details of these effects depend greatly on the character of the 
speed of sound profile: If the sound does not increase sufficiently rapidly 
above and below the axis, for example, the time stretching may not occur, or 
may even be reversed, with the sound traveling straight down the axis arriving 
first. Additionally, the exact location of the source and receiver with respect 
to the channel axis is very important with respect to the envelope of the 
received signal and the amount of time stretching encountered. For the 
typical profiles in the northern hemisphere, the increased duration of the 
received signal over that sent can be estimated by 



if source and receiver are on the channel axis, and less if they are not. 
The quantity t is the total time of flight between source and receiver and 
Ac is the maximum variation in the speed of sound between the upper and 
lower boundaries of the channel indicated in the previous page. Typically, 
this leads to about 9 seconds of stretching for each 1000 nautical miles 
the signal has travelled. 



At . 1 Ac 

t 3 c 
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NORMAL MODES 



Another approach to the propagation of sound in the ocean is 
Normal Mode Theory. If the mathematical difficulties of this 
approach can be overcome, it yields more exact values for the 
transmission loss than can be obtained using Ray Theory. 

Normal Mode Theory yields a sum of individual solutions 
(each called a normal mode) which must be combined to obtain 
a complete expression for the sound field. 



These special factors Z n are called "normalized eigenfunctions" 
They are found by solving the wave equation and specifying the 
behavior of the pressure at the ocean surface and bottom. The 
solution for each Z n also gives a required value for c n . Both 
the depth dependence of Z n and the value of c n depend on the 
frequency of the sound waves. 

This kind of a solution can be obtained whenever the speed of 
sound profile c(z) is shaped so that the sound can be trapped 
in a channel. The trapping of sound in the mixed layer is one 
such example, and the propagation of sound to long ranges in 
the SOFAR channel is another. 




Each one of these normal modes has its own behavior 




This part looks just like a cylindrical wave 
traveling outward with its amplitude falling 
off like l//r* . Notice however that it has 
its own special speed c n . 



This part is an amplitude factor which changes with 
depth. Notice that each normal mode has its own 
special depth-dependent amplitude factor. 
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NORMAL MODES FOR A CHANNEL 





The above sketches represent the normalized eigenfunctions of the 
first three and sixth normal modes for sound trapped in the 
SOFAR channel. These sketches show that the greatest part of 
the pressure field is found in the channel, but the long "tails" 
which extent to depths far below the channel axis show that 
there is sound ensonifi cation at depths greater than those 
predicted by drawing ray paths. 
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CUTOFF FREQUENCY AND LOSSES 







i 

i 



I 





Each normal mode has its own cutoff frequency f . If the sound 
has a frequency below that of cutoff for a particular mode, then 
that mode cannot be excited to carry sound energy in the channel. 

If the frequency of the sound is slowly increased with time, the 
above sketches show how a normal mode can be excited. Just above 
cutoff, the normal mode has a very long tail which goes to great 
depths beyond the edges of the channel. As the frequency increases 
further above cutoff, the tail gets shorter and shorter, until 
at frequencies well above f there is very little sound existing • 
outside the channel. 

This means that for a specified frequency, sound can be trapped and 
sent to large distances only in those normal modes whose cutoff 
frequencies are less than the frequency of the sound. Normal modes 
with higher cutoff frequencies will have longer tails than those 
with lower cutoff frequencies, and those whose cutoff frequencies 
exceed the frequency of the sound cannot be excited at all. 

The larger the tail extending out of the channel, the more energy 
can be lost from the sound in the channel. This means that 
higher modes (those with higher cutoff frequencies) will attenuate 
faster and be less significent at great distances. This effect is 
sometimes called mode stripping . 

These energy losses for each mode can be handled by specifying an 
absorption coefficient for each normal mode. 

where oc n increases with increasing n. 



Then, we have 

e 
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TRANSMISSION LOSS (NORMAL MODES) 



n 





If incoherence 

P= + ■■■ 




low frequency- 
short ranges 
homogeneous water 
smooth boundaries 



high frequency 
long ranges 
inhomogeneous water 
rough boundaries 
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Transmission Loss in dB 



RAYS, NORMAL MODES, AND REALITY 

Normal mode theory is a complicated mathematical formalism which does not 
lend itself to simple intuitive interpretation. However, it is an exact 
theory which avoids certain approximations inherent in ray theory. It is 
useful in discribing situations for which ray theory fails. An example 
of one such situation is illustrated below. 



Range 




The source is in the layer and the receiver is below the layer. Ray 
theory predicts a shadow zone at large enough ranges. While there can 
be some scattering from the ocean surface into the shadow zone, this 
will be small for low frequencies and low sea states. The surface 
interference pattern at short ranges shows that the sea was quite 
smooth (the observed sea state was about 1 or 2). 




Notice that the experimental data show significent ensonification of the 
shadow zone. This cannot be predicted by the ray theory, and the scattering 
of sound from the surface is not strong enough to account for it. 

Reference: Pederson and Gordon, Jour. Acoust. Soc. Am. 37 , 105 (1965) 
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SUMMARY OF TRANSMISSION LOSS EXAMPLES 



For both source and receiver in the surface layer, and at ranges 
greater than the transition range, the transmission loss can be 
estimated by 



TL • 10 log r. + 10 log r + ar + br/r . 

t s 



For a convergence zone detection, we can write 
TL = 20 log r + ar - G . 



If the propagation path is bottom bounce, and the angle of depression 
of the sound beam in that mode is sufficiently large, then the simple 
geometrical model gives 



TL 



20 log 



r 

cos 0 



+ 



ar 

cos 9 



+ 



BL . 



For propagation between a deep point and a shallow point over the 
reliable acoustical path, an approximate transmission loss formula is 



TL 



20 log 



ar 



cos 0 



cos 9 



Finally, for the propagation of sound trapped in the S0FAR channel, 
if the range is greater than the transition range r^.’, the transmission 
loss can be estimated by 



TL 



10 log r + 10 log r^' 



+ ar . 
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COMB ININS SIGNALS 



Now, let us no longer restrict ourselves to tonals (monofrequency 
signals). Because we really don't deal with only a single frequency, 
we have to find out how signals of different frequencies combine. 
Consider the combination of a large number of monofrequency signals j 



\/^ 



The combination of many signals, each a pure tonal, will result in 
a total signal which does not appear to have definite amplitude 
or period. However, the composite signal is an acoustical wave, 
and the total acoustical energy it possesses is the sum of the 
acoustical energy of each of the individual tonals. Thus, 
if E>i> is the total energy of the composite signal, we have 

N 

Erp = + E2 + * * ' = £ E^ 

i = I 

where E^ is the energy of signal 1, E 2 that of signal 2, and so on. 

For the complicated signals encountered in real life, which have 
large variations from moment to moment in their amplitude and 
shape, the frequencies of the individual tonals which combine to 
form them are extremely closely spaced. Indeed, it is often 
more convenient to talk about the energy contained in the 
collection of tonals which lie between some specific frequency and 
another frequency which is only inf initisimally different. 

If aE is the total energy of all the tonals lying within the 
incremental bandwidth aw, then we can define the quantity £ 
as the average energy density over the incremental bandwidth aw by 




A MS 



Notice that the energy density has the units of energy/Hz. If we 
analyze some complicated signal, therefore, we can describe it 
in terms of the energy density it contains in each incremental 
interval of frequency over the entire range of frequencies which 
comprise it. 




N 



• ■ 

vvwvvwW* j 
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ENERGY SPECTRUM AND PRESSURE SPECTRUM LEVEL 

An example of the analysis of some complex signal into the 
energy densities found in each frequency interval is suggested 
in the following curve i 




The total energy of the signal is given by 

E T ~ £ Amt + £ 2 ■+ * * • - ^ £ Jur i 

i 

which is the area under the curve of energy density vs frequency. 

If all of the ^w^ are chosen to be 1 Hz wide, then the graph gives 
the energy to be found in each 1 Hz interval, and the curve is 
called the energy spectrum of the signal. 

It is more conventional in acoustics to describe the signal not 
in terms of the energy spectrum, but rather in terms of the 
pressure spectrum level PSL, which we shall define here as 

where the quantity p , the effective pressure per Hz, generates 
the associated energy found in the 1 Hz bandwidth. For example, 
if we receive a certain acoustical signal, pass it through a 
filter with a 1 Hz bandwidth, and obtain the effective pressure 
amplitude of the resultant filtered signal, say 100 ^.Pa, then the 
value of p for that 1 Hz bandwidth is 100 ^ Pa/l Hz * 100 ^-Pa/Hz. 

A representative curve of the pressure spectrum level vs frequency 
is presented below s 
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AVERAGED ENERGY SPECTRUM AND AVERAGED PSL 



Very often, what is presented as the energy spectrum or the 
pressure spectrum level is not that which would be measured with 
a 1 Hz filter, but is rather calculated from measurements made 
with a filter of larger bandwidth. 

The overall energy to be found in some arbitrary bandwidth 
w is the area under the energy spectrum vs frequency curve 
over the appropriate range of frequencies; 




This gives the total energy E w found within this bandwidth, and 
does not say how e is distributed over w. We can therefore write 
the averaged energy spectrum £ as £ = E w /w, A plot of £ vs 
frequency would give a smoothed curve, lacking therdetail of the 
true curve, but sufficient to calculate the total energy present 
over bandwidths comparable to or larger than the bandwidths selected 
to obtain £" . 

Analogously, it often happens that what is presented as a 
pressure spectrum level is in reality a smoothed curve, the result 
of plotting data obtained with filters having bandwidths appreciably 
larger than 1 Hz t 




Use of too broad a bandwidth in calculating averaged spectrum 
levels will wash out details. The spikes in the true spectrum 
level in the above example are lost to view because the bandwidth 
over which the averaged levels are taken is so wide that the energy 
from the spike is small compared to all the energy from the 
background in each bandwidth. 
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THREE WAYS OF DOING IT 




r irvr . 

Spectrum Leve'|^f 
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The first graph presents the true 
spectrum level and the other three 
present the same spectrum, except 
that the spectrum levels have 

been averaged over the bandwidth 
w. Graph 2 presents the 
averaged spectrum level in a 
bar graph format. Graphs 
3 and 4 connect the data 
with either straight-line 
segments or a smooth 
curve. Notice that 
Graph 4 gives no 
information relating to 
the bandwidth of the filter 
used to obtain the averaged spectrum 
levels? this would have to be 

indicated separately, and is not 
likely to be done. Graph 3 is 
somewhat better, but the use 
of the straight-line segments 
seduces the eye into assuming 
that the true spectrum level 
! closely follows the 

i straight-line segments. 

Graph 2 is the least 
subject to being 
misinterpreted: The 

width of the individual 
bars gives the bandwidth 
of the filter quite 
visually, and each bar strongly 

emphasises by its flat top that the 
value of the PSL appropriate to 
that frequency interval is the 
averaged value. 

Because of the fact that 
the averaged spectrum level 
smooths out and erases 
detail which would 
otherwise be found in 
the true spectrum level, 
it is very important 
that the extent of averaging 
over frequency be well- 



Frequency 






Frequency 




indicated. This is accomplished quite 
easily and unambiguously by the format 
of either Graph 2 or Graph 3. Further, 
Graph 2 clearly indicates the 
averaging of data by the flat-topped 
nature of the curve, whereas 
Graph 3 can tend to mislead 

the eye in a casual examination. 
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SPECTRUM LEVEL AND BAND LEVEL 



Recall that the total energy in a bandwidth w was the area 
under the energy density vs frequency curve. The analogous 
quantity in terms of pressure levels is the band level BL. This 
is the expression in dB terms of the total effective pressure 
generating the acoustical energy in the specific bandwidth. 

NOTICE THAT THE BAND LEVEL IS NOT THE AREA UNDER THE 
SPL vs FREQUENCY CURVE. 



This is because the log of a product is the sum of the logs of 
the individual factors. 

Consider a plot of energy density £ vs frequency as shown* 



£ , 

£ 

o 


k 




Mr — v 


^ 







-f 

It is clear that the total energy E in the bandwidth w is given by 

E = £ o w 

If we express this in log form, we have 








-4- 



Q 



It turns out in acoustics that taking 10* log(E/E re; j.) is the same 
as taking 20* log(P/P re f ) . This means that the above equation 

is the same as 



loA 



_P_ 



* lo JZ>} £ir— + 10 £ 
7 V t ,t 



1 



1 W ■ 



or, treating the quantity 1 Hz as understood but suppressed, 



BL = PSL +10 log w 

In general this expression is exact only if the spectrum 
level is constant over the bandwidth w. However, it is 
reasonably accurate if the true spectrum level does not 
fluctuate by more than a couple of dB over the range of w. 



Notice that this formula can be turned around* If a band 
level BL has been experimentally obtained by measuring the 
acoustical energy over some bandwidth w, then the averaged 
spectrum level over this bandwidth can be found from 

PSL(averaged ) = BL - 10 log w 
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COMPARISON OF A SPECTRUM LEVEL AND A BAND LEVEL CURVE 
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Because of the (10 log w) term which must be added to the 
PSL.to obtain the appropriate BL for a certain frequency interval, 
it is seen that the shape of the BL curve is such that the BL 
is stronger at higher frequencies. This is because the BL measures 
the total acoustical energy in some frequency interval, whereas the 
PSL measures the energy per unit Hz. 

Both kinds of curves have their advantages and disadvantages. 

A curve of the true spectrum level gives the maximum possible 
amount of information, but requires considerable mathematical 
analysis if it is desired to obtain the band level within some 
specified frequency interval from the spectrum level. A curve 
of averaged spectrum level is not quite so difficult (timeconsuming) 
to use, and band levels from it can be found more easily, but often 
important details are lost. A curve of band level is quite direct 
to use, as long as the bandwidth agrees with the particular receiving 
system being used, but it also has smoothed out much detail. 
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A SAMPLE CALCULATION 



As a simple example, let us obtain the overall band level 
for the PSL curve presented below: 



From the steplike nature of the curve, we can specify two bandwidths 
w^ and wg as shown, with the associated averaged spectrum levels. 

For each of these, we can calculate the appropriate band level: 



So far, we have had to do no more than a little curve-fitting and 
application of formulas which have been presented previously. 

Now, however, a new problem has appeared. We must combine the 
two band levels into a single bandlevel for the total frequency 
interval w^ + W£. The combined BL is not the sum of the individual 
ones, because we are in dB: Energies add but logs of energies 

do not. What must be done is to determine the amount of energy 
given by BL(w^ ) and the amount of energy given by BL(wo), and 
these energies added together and then expressed as a band level. 

Let us illustrate this with an example. We shall first do 
the example in the method indicated above, and then present a 
simple nomogram which avoids all the work. 
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BL(w^ ) = B(w^ ) + 10 log w^ 
BL(w 2 ) = B(w 2 ) + 10 log W 2 
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NOMOGRAM FOR COMBINING LEVELS 




As an example of the application of this nomogram, recall the 
problem worked out on the previous page. The two band levels 
which were calculated there were BL^ = 170 dB re l/*Pa and 
BLp = 168 dB re l^Pa. The difference is +2 dB? finding this 
value, projecting upward, and reading to the left gives 2.1 dB, 
so that the combination of these two levels yields a total 
band level of 170 + 2.1 = 172.1 dB re l^Pa. Since the usual 
uncertainty in levels exceeds a few tenths of a dB, it is 
appropriate to round this combined level to 172 dB re l^fa. 



TONALS COMBINED WITH BACKGROUND 




-1 

Vi 
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Frequency 




When there are strong tonals present, the band level cannot 
be expressed as B + 10 log w, because of the additional energy 
contributed by the spike(s). We must calculate the energy from 
the background spectrum (without the spikes) and add this to the 
energy of the monofrequency spikes which are also present: 

BLg * B + 10 log(w-l ) = B + 10 log w if w^l 

BL A = A + 10 log 1 = A. 

The total band level is found by combining 3L A and BL g with the 
help of the nomogram. 

Notice that if the tonal is relatively strong and the bandwidth 
narrow enough, the combined level will be due essentially to the 
tonal. As the bandwidth w increases, however, the influence of 
the (10 log w) term will become stronger, causing BLg to increase. 
Eventually the contribution from the background will get sufficiently 
larce that BLg will be sufficiently large that the influence of 
BL a will be lost. 
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THE SONAR ’SQUAT! ON 



EL ^ML + DT 

i 



EL 55 Echo Level for active sonar or the received Signal Level for passive 
sonar. It is the level in dB of the desired signal vhen it reaches 
the receiver. 

11L = Masking Level . It is the level, in dB of the undesired signal at the 
receiver which is competing vith the Echo Levol. 

DT * Detection Threshold . This is the number of dE by vhich the Echo Level 
must exceed the Masking Level to allov the receiving system to 
register a detection vith a specified degree of confidence. This 
quantity involves the probability of detecting the target in the 
presence of noise and also the probability of giving a false alarm. 



The SCIiAFi equation states that signal at the receiver must equal or exceed 
background noise by an amount specified by the Detection Threshold if a 
specified degree of performance is to be achieved. 
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LIST CF SYMBOLS 



EL 


■ 


Echo Level 


ML 


m 


Masking Level 


DT 


m 


Detection Threshold 


SL 


m 


Source Level 


TL 


m 


Transmission Loss 


PSL 


m 


Pressure Spectrum Level 


BL 


m 


Band Level 


SSL 


m 


Source Spectrum Level 


AG 


m 


Array Gain 


NL 


m 


Noise Level 


NSL 


m 


Noise Spectrum Level 


DI 


m 


Directivity Index 



w * bandwidth 

d - detection index 

wt * time -bandwidth product 
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PASSIVE SONAR 




SL - TL ML + DT 



where SL * Source Level of the target . 
TL = Gne-vay transmission loss 



dach term in the above equation will now be discussed separately. 
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(sl)- TL £ ML + DT 
SOURCE LEVEL - PASSIVE SOLAR 



T . Tonals (monfrecuency components) 

Sources: Propulsion machinery (Main motors, reduction (.ear) 

Auxiliary machinery generators, pumps, air-conditioners) 
Propeller-induced resonant hull excitation 
Lydrodynamically-iriduced resonant excitation of cavities, 
plates, and appendages 



SL = Level of acoustic energy 1 meter from tie acoustic center of the 
source. 

measurements are usually made at a distance from the source and 
extrapolated back to 1 meter. 



II . broadband noise 

Sources: Cavitation at or near the propeller 

Radiated flow noise 
Caviation at struts and appendages 



SL * SSL + 10 log v 

where w is the bandwidth of the system and 

SSL = pressure spectrum level 1 meter from the acoustic center 
of the source. 

If the SSL is not essentially constant over w, the bandwidth 
must be subdivided into smaller segments, the LL of each 
found, and then all the LL's combined to give SL 
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TYPICAL SOURCE SPECTRUM LEVELS 



Typical source spectrum levels for a conventional submarine on electric 
drive at periscope depth. SSL re l/>Pa. 




Adapted from Urick, Principles of Underwater Sound for Engineers 



SL - TL DT 

MASKING IEVEL 







Unwanted noise may come in uniformity from all direction or it may arrive 
from a prefered direction. 



ML - NL - AG 

where 

ML ■ masking ictel 

NL - noise level in the environment 

AG ■ array gain, the improvement (reduction) in dB against the 
noise field resulting from the array design. 
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where 



and 



Note: 



ML - ML - AG 
Koise Level 



ML = NSL + 10 log w 



NL = noise level, the acoustic energy in dB found in the bandwidth 

of the system , 

NSL * noise spectrum level, pressure spectrum level of the unwanted 

noi se , 



w * the handwidth of the system in Hz. 



The same limitations apply in using this formula as were 
pointed out in the discussion of the broadband SL. 
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NL * NSIJ + 10 log v 



NOISE SPECTRUM LEVEL (Deep water) 




Adapted from: Venz, Jour. Acoust. Soc. Am. 3U , 1936-1956 (1962) 

Perrone, ibid 1*6, 762-770 (19697 
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N0IC3 



KL = (IJSL + 10 lot vj 
: SPjXTKUM LEVEL (Self -noise) 



Sources of self-noise: Hydrodynamic noise due to flov past receiver 

Acceleration of receiver 
Machinery noise of ovn "ship" 



Values usually presented as an "equivalent omnidirectional" noise spectrum 
level. 



KL( self -noise) * HSL( self -noise) + 10 lot w 



If .significant self -noise is present, it must be combined vdth the ambient 
noise level to obtain the total noise level. Since these two sources of 
noise are independent, their energies add and the nomo£ram can be used. 
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ML 



NL - f AG 



ARRAY GAIN AND DIRECTIVITY INDEX 



Array gain is the decrease in masking level from the noise 
level brought about by the design of the receiving array* 



If the noise is isotropic (arriving uniformly from all directions) 
and the signal arrives from one direction only, then 

AG = DI 



where 



DI = Directivity Index 

= 10 log (the ratio of the noise energy sensed 

by an omnidirectional receiver to that 
sensed by the directional receiver) 



In this case, 



ML = NL - DI 



Array gain can include both the design of the receiving array 
and also the signal processing. Signal processing such as correlation, 
by taking into account the different spatial or temporal properties 
of the signal and noise can result in an array gain different from 
the directivity index. 
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Pair by pair, 
the signals 
cancel 




axial' 



DI - 0 dB 



Omnidirectional receiver 



DI • 20 dB Receiving only 1$ of the energy present in the noise 
field. 

1 depends on site of the array compared to the wavelength. For example, 
or a long line array of length L much greater than ^ 

DI - 10 log (2L/<\ ) 

morally, as frequency increases the receiver becomes more directive (DI increases) 
at, if the frequency is too high we may get more major lobes. 
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A PHYSICAL MEANING OF DIRECTIVITY INDEX 




Surface area of the unit sphere * IX (>) 

Area of the illuminated patch * 

where JX is the solid angle subtended by the patch. 

The fractional area of the sphere which is illuminated 
is given by SL /4tT and, for this simple example, the directivity 
D is given by n 

The directivity index is the expression of the directivity D in 
dB terms, 



DI = 10 D 



A more general definition of the directivity is given as 
follows: If the pressure field generated by an array acting like 

a source (or the relative sensitivity of an array acting like a 
receiver) is written in the form 

P(r,0 ,0) - P ax (r) H(e,0) 

where P ax is the pressure on the acoustical axis of the array, 
then H(9,f0 is normalized to have a maximum value of unity on 
the acoustical axis, and to be no greater than unity for any other 
direction. The directivity can then be expressed as 
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SOME ROUGH ESTIMATIONS OF DI 



If we have a line array with a single major lobe, and if 
the angular width of the major lobe is 9 (in radians), then if 
the array is surrounded by a unit sphere, the illuminated patch 
appears as below 




and the directivity can be crudely estimated to be 

±311)1 Z 



2irO)-9-(i) 



e 



( 0 in radians ) 



where it must be understood that 9 is small. 



A little more careful estimation of the directivity from 
the geometry of the radiation (or reception) pattern leads to 
some relatively simple but only approximate estimates for the 
directivities of arrays of simple design, 

Line- like arrays? 



r\ 7[ Ltiqth 

u l A 



» I 



Piston-like arrays: 








»l 



These formulas can be used as rough estimates if the arrays are 
highly directive (DI in excess of about 6 dB at the least), are 
designed to have only one major lobe, and are not strongly shaded 
or steered. 

An array is shaded when different elements have different amplifications. 
This is done to alter the width of the major lobe and the relative strengths 
of the side lobes. Narrowing the major lobe will increase the strengths of 
the side lobes, and vice versa. An array is steered to point the major lobe 
in various directions with respect to the array. This is done by introducing 
time delays in the outputs of the elements. 
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PASSIVE SONAR EQUATIONS 



Lat's put it all togathar: 

BROADBAND WITH NEGLIGIBLE TONAIS 

SSL + 10 log v - TL > NSL * 10 log w - AG ♦ DT 

SSL - TL ^ NSL - AG ♦ DT 



TONALS WITH NEGLIGIBLE CONTRIBUTION FROM BACKGROUND 



SL - TL > NSL + 10 log v - AG ♦ DT 
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FIGURE OF MERIT 



(FOM) 



The passive SONAR equations can be rewritten by solving forTL: 



TL ^ POM 

where 

POM - SL - (NSL + 10 log w - AG + DT) 
or 

POM - SSL - (NSL - AG + DT) 
for total or broadband detection respectively. 
Graphically: 





ranges for which detections are possible at conditions 
equal to or better- than specified by DT. 
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To take advantage of a narrow bandwidth and still permit reception of a 
Doppler shifted signal, it is necessary to design a system that employs 
many parallel filters which, while individually narrow, cover the entire 
frequency range of interest. 

The output of such a system might be presented in this manner: 




At any instant, the output of each of the individual filters might look 
like this: 
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Probability Density of Amplitude A 



SL - TL < ML + @ 
DETECTION THRESHOLD 



The detection process consists of designating a threshold which, when exceed, 
causes a detection to be recorded. If the signal is much stronger than the 
noise, it is clear that a threshold can be defined that will allow valid 
signals to be recorded while r ignoring the noise. However, when the 

signal and noise are of comparable size, any threshold that will catch a 
reasonable number of valid signals will also record "detections" when a valid 
signal is absent. 

The probablity of obtaining a certain amplitude of the noise alone as well as 
for the signal- and- noise is sketched below: 




where 

A n is the mean value of the amplitude with noise alone 
A sSn is the mean value of the amplitude with signal and noise 
cr is a measure of the width of the curves (assumed equal) 

T is the threshold 

The shaded area marked "detection" is the probablity of making a valid detection 
when the threshold is T. (The area under the curve is taken as unity.) And 
the shaded area marked "false alarm" is the probablity of making a "detection 
when the signal is absent. 

Let P(D) ■ the probability of detection 

P(FA) ■ the probability of a false alarm 
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DETECTION INDEX 







Amplitude 




We need to find some simple 
criterion which will describe 
how the curves of noise alone and 
signal and noise combined overlap, 
since for a given value of 
threshold, it is the overlap of 
the curves which establishes 
the relative values of P(D) and 
P (FA) . 

Notice the configuration of the 
two curves given in the first 
graph. The curves are quite well- 
separated, in that their peaks 
have amplitudes whose difference 
is larger than the o' . These 
two curves could be considered 
well-eeparated . 

In the second graph, each curve 
has the same average value of 
A as before, but now the <T * is 
twice as large. The curves now 
are not well-separated. 

In the third graph, the curves 
have the same O'' as in the second, 
but the two have been shifted 
further apart so that they are 
well-separated. Indeed, compariso: 
of the first and third graphs 
will reveal that the amount of 
overlap of the pairs of curves 
is the same. Notice that the 
ratio of the difference in the 
verage values of the amplitudes 



— t 1 *- 

A v A 




Aa&n 3114 t0 the ^ tile 

first graph is identical with 
that for the third. Thus, it 
is the separation between the 
two curves ih units of the width 
of the curves which measures 
how well-separated the two curves 
are. It is thus plausible to 
define the detection index d by 

- ■ 
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Let us study the effects of changing the positions of the threshold and the 
detection index on the probabilities of detection and false alarm. If the 
original configuration of the noise curve, signal and noise curve, and 
threshold are as given in the top curve, then we can establish the indicated 
P(D) and P(FA) by the shaded areas. Now, if the threshold is lowered, as in 
the second figure, then P(D) shows some increase, but the P(FA) increases 
rather drastically: lowering the threshold will increase the probability of 
a detection, but the "clutter" will also increase. On the other hand, if the 
noise curve and threshold are maintained constant, but the signal gets stronger, 
then the curve of signal and noise will move to the right. This means that the 
detection index increases, and the P(D) increases while P(FA) remains the same. 
This says no more than that a stronger signal has a greater probability of 
being detected. 

If the statistical properties of the noise and the combined signal and noise 
are known, then it is possible to make curves of the P(D) and P(FA) relationships 
for various values of detection index. These are known as Receiver;- Operating 
Characteristics (ROC) curves. 
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P(D) 



RECEIVER OPERATING CHARACTERISTICS 
ROC 




P(FA) 

Typical ROC Curves 
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FALSE ALARM RATE 



The rate of false alarms is the product of the sampling rate and P(FA). 
Since sampling rates are generally high, the importance of a very low 
P(FA) should be noted. 

For examples If the output of each of 10 parallel filters is sampled 
for one second each, a P(FA) of 1$ will result in an average of one 
false alarm every 10 seconds. 

The sketch below shows a. representative plot of the relationship between 
d and P(FA), for a P(D) - 0.5*. 




P(FA) 
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ENERGY DETECTION 



For a passive detection system, 

a completely unknown signal, and Gaussian noise, best performance is 
obtained by filtering in a bandwidth w, taking the square of the signal, 
and accumulating over a processing time t. If the signal-to-moise 
ratio is small, it can be shown that the detection threshold for this 
process is 



DT e • 5 lo£ (d/wt) 



where wt is called the time-bandwidth product. 

To decrease DT for a fixed d, it is necessary to increase the time -bandwidth 
product. 



Footnote: 

The detection threshold used in this material is that calculated over 
the complete bandwidth of the processor. It often occurs that what is 
specified is the detection threshold in a one Hertz bandwidth, DT(l). 

This latter quantity occasionally leads to SONAR equations which have a 
slightly simpler format when it is assumed that the spectrum levels are 
constant over the bandwidth of the system. Otherwise, we feel that the 
use of nr(l) introduces more possibilities for confusion, so that its use 
creates more problems than it solves. (In particular, if Urick's excellent 
book is studied, it must be fully understood that whereas he uses DT throughout 
most of the text, his chapter devoted to the detection threshold uses DT(l).) 
The formula relating the two quantities is 
« 

DT(1) - DT + 10 log w 
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DOPPLER SHIFT 



C=3 




AJ -- 2 %$ 



While decreasing w will improve the ability of a paasive SONAR to detect 
tonals, it cannot be made smaller than the expected Doppler shifts. 

For a passive SONAR system, this requires that the minimum acceptable 
bandwidth be 

w(in Hz) = 0.7 x (frequency in kHz) x (range rate in knots) 
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EXAMPLE 1 



A submarine with a 160 dB re l/<Pa line at 500 Hz crosses a convergence 
zone at 1*0 km from an omnidirectional hydrophone. The bandwidth of the 
receiver is 100 Hz and the sea state is 3« Assuming a convergence gain 
of 12 dB, find out how long the signal must be accumulated to give a 
P(D) - 0.50 and a P(FA) - 10“^. 



Solution: 

TL = 20 lot r + ar - G = 20 log (1* x 10^) ♦ 5 x 10 ^ • U X 10^ - 12 * 82 dE 

ML » NSL + 10 log v - DI - 66+20-0 « 86 dB 
SL - TL > ML + DT 

DT ^ SL - TL - ML = 160 - 82 - 86 - -8 dB 

DT ■ 5 log (d/vrt) 
where d • 15 

log (vt/d) ■ 1.6 
wt/d • itO 
wt • 600 
t - 6 sec 



88 



EXAMPLE 2 



Sams as Example 1, but let w be Just vide enough to encompass a 20 knt 
target. 

Solution: 



where 



V 


- 


0.7 x (frequency in kHz) x (range rate in 


knots) 




- 


0.7 h 20 - .. 7 Hz 




ML 


9 


NSL + 10 log v - DI - 66 + 8.5 - 


74.5 dB 


DT 


< 


SL - TL - ML - 160 - 82 - 74.5 - 


3.5 dB 


DT 


- 


5 log (d/wt) 




d 


m 


15 





log (d/vt) * 0.7 



d/vt » 5 

wt - 3 

t ■ 3/7 sec 
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EXAMPLE 3 



A conventional submarine at periscope depth is traveling at 4 kts. The 
Sea State is 3 and the layer depth is 100 m. At what range will the 
submarine be detected by a 36 m deep hydrophone listening at 1000 Hz 
if DI = 20 dB and DT = 0 dB? 

Solution: 

SSL - TL ^ NSL - DI + DT 

TL < SSL - NSL + DI - DT = 120 - 62 + 20 = 78 dB 

For sound to be trapped in the duct, 

f > f . = 2-10 5 /D 3/2 = 200 Hz 

r mm 

So we use the TL for a mixed layer 



TL = 


10 log r + 


10 log r + ar + b(r/r ) 

s 


ri- 

ll 


1.05-10 2 D/(d 


1/2 2 

- z ) ' = 1.05-10 (100/8) = 1310 m 

s 


r = 


8.4-10 2 D 1/2 


= 8.4* 10 2 * 10 = 8400 m 



s 

From the graphs and equations 
a = 6*10 5 dB/m 

b = 1.7 dB /bounce 

Notice that ar + b(r/r ) + (a + b/r ) r 

s s 

Thus 78 = 10 log 1310 + 10 log r + 

10 log r + 2.6-10 -4 r = 47 

Now, trial and error 



r 


10 log r 


2.6*10~ 4 r 


TL 


o 

i — i 

in 


47 


13 


60 


o 

i — i 

CM 


43 


5.2 


48 


4 

1.5-10 


41.8 


3.9 


46 


4 

1.7-10 


42.3 


4.4 

1 


47 



r = 1.7-10 4 m 



= 2.6-10 -4 r 

2.6-10 -4 r 
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EXAMPLE k 



Repeat Example 3 but with DI * 0 dB. 

Solution: 

TL £ SSL - NSL + DI - DT • 120 - 62 • 58 dB 

-4 

TL » 10 log r + 10 log r + 2.6 x 10 r 

t 

58 - 10 log 1310 + 10 log r + 2.6 x lO -4 r 



58 - 31 ■* 


10 log r + 


2.6 x 10“ 4 


r 


10 log r + 


2.6 x 10 -4 r 


= 27 




r 


10 log r 


2.6 X 10“ 4 


TL 


5 x 10 2 


27 


0.13 


27 



But r < 

So we must use this form of the TL.' 

TL * 20 log r ♦ ar 

58 - 20 log r 

log r - 2.9 

r ■ 800 m 



Bingo 
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EXAMPLE $ 



Repeat Example 3 with f » 100 Hz. Assume "heavy shipping" as would be 
encountered in the Mediterranean Sea. 



Solution: 



f < f , 



min 



No trapping in the layer. Use 

TL » 20 log r + ar 

TL £ SSL - NSL + DI - DT ■ 133 - 76 + 20 - 77 dB 
20 log r + 10** 6 r - 77 



r 


20 log r 


10" 6 r 


TL 


8x 10 3 


78 


8 x 10“3 


78 


7 x 10 3 


77 


7 x 10’ 3 


77 



r * 7 x 10 3 m 



Bingo 
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ACTIVE SONAR 




SL - 2TL + TS \ ML + DT 

\ 



where 

SL ■ source level of transmitter 
2TL ■ two-way transmission loss 
TS ■ target strength 



Each term in the above equation will now be discussed separately. 



Note; the above active SONAR equation is for a monostatic case, 
wherein the source and receiver are coincident so that the TL 
from source to target is the same as the TL from target to 
receiver. Our discussion will usually be restricted to this 
case. Extension to the bistatic situation in which source and 
receiver are at different locations is trivial: If the 

transmission loss from source to target is TL, and from target 
to receiver TL' , then all that is required is to replace 
2TL with TL + TL' , so that we have 



SL - TL - TL’ + TS > ML + DT 
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LIST OF ADDITIONAL SYMBOLS FOR ACTIVE SONAR 



SL = source level of the active sonar system 
TS » target strength of the desired target 
RL * reverberation level 

TS' * target strength of the scatterers causing reverberation 
TL = transmission loss from source to target 

TL' = transmission loss from target to receiver, equal to TL in 
the case of monostatic systems 



L e = "effective length" of a target 

t * pulse duration of the acoustical signal from the source 
s y = a measure of the scattering cross-section per unit volume 
s g ® a measure of the scattering cross-section per unit area 
S v 35 scattering strength for volume scatterers 

S_ = scattering strength for surface scatterers- 

s 
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@) - 2TL + TS ^ ML + DT 
SOURCE LEVEL - ACTIVE SONAR 



SL • sound pressure level 1 m from the acoustic center of the 
transmitting transducer. 



SL ■ 1© log (Power) + DI + 171 



where 



or 



where 



and 



Power • Acoustic power output of the transducer in Watts 

DI = 10 log (the ratio of the intensity on the 

acoustic axis to the intensity at the same 
distance from an omnidirectional transmitter 
with the same acoustic power output) . 



DI 

P 0 (r) 



20 log P(r)/P Q (r) 

pressure produced by an omnidirectional transducer 
with the same power output 



171 ■ factor to take care of the units 



4 

*4 


2i|0 


•V 




H 




2 


220 


€ 




fH 

$ 


200 


© 

►4 




© 

o 


180 


i 




CO 


160 




Acoustic power output (Watts) 
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SL - 2TL + (TC>)^ ML + DT 
TARGET STRENGTH 



TS • (sound pressure level of reflected wave 1 m from the 
acoustic cenfefcr of the target in a direction back 
towards the receiver) - (sound pressure level of 
the incident wave at the position of the target if 
the target were absent). 

or 

TS - 20 log P^r'-lJ/Ptfr) 

where 

P'(r'*l) is the pressure of the reflected wave 1 m from the 
acoustic center of the target in a direction back towards the 
receiver. 



Target Strength depends on shape, construction, and orientation of the 
target, and on the frequency and pulse length of the sound wave. 
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EFFECT OF TARGET SIZE 




Because of the effective length L of the target, portions of the 
incident pulse which reflect froiri the tail (in the above sketches) 
travel a total distance of 2L further than the reflections from 
the nose. This means that thi apparent duration of the pulse 
is stretched from t to (t + 2L /c) . For very long pulses, 

the additional time interval corresponding to the round trip, 

2L /c, can be neglected. For short pulses, however, this stretching 
can be quite important in determining the target strength. 
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The Table shows a few representative, but very approximate, values of 
target strength for a variety of objects. 



Target 


Frequency 




in kHz 


Fleet Submarine 


10 


(approx 100 m long) 


2h 


S Class Submarine 
(approx 70 m long) 


2k 


Torpedo 


2k 


Surface v essel 




3 ft mine 


30-90 



Approximate Target Strength (dB) 
Beam Aspect Bow-Stern Aspect 

35 30 

25 16 

20 13 

-20 

10 to 30 
-8 



Extracted from: Physics of Sound in the Sea, NDRC 

Principles of Underwater Sound for Engineers, Urick 
Principle and Applications of Underwater Sound, NDRC 
NAVSHIPS 0967-129-3010 
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SL 



2TL + TS (ML) + DT 



MASKING LEVEL 

The undesired noise competing with the desired signal can 
arise from two possible mechanisms in the case of active SONAR: 

Masking from ambient or self noise 

Masking from reverberation 

Which of these two mechanisms is the dominant depends on the 
range between source and target. Generally speaking, for low- 
powered active systems the performance will be noise-limited, 
in that the maximum effective detection range will be determined 
when the echo level diminishes below that level for which it can 
be extracted from the ambient or self noise. For higher-powered 
active systems, however, reverberation becomes an important source 
of interference. Reverberation tends to decrease with increasing 
range, but less rapidly than the echo level. Thus, if the echo 
level has diminished until it is buried in the reverberation, 
but the reverberation has not yet fallen sufficiently below 
the ambient (or self) noise, then the maximum detection rangp 
is determined by the reverberation. The two situations are 
suggested in the figures: 




The choice is a little more subtle than suggested by the simple 
curves above. The detection threshold for noise masking may 
be different than that for reverberation masking, And the 
detection threshold for reverberation-limited performance may 
depend strongly on the doppler- induced difference between the 
frequency of the reverberant sound and that reflected from the 
desired target. 
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AMBIENT AND SELF NOISE 



The masking level is given by 

ML = NL - AG 

where the noise level NL is that from the ambient noise of the 
ocean, noise generated by the receiving platform, or the combination 
of the two if they are of nearly equal strengths. 

Recall that if the noise spectrum level is sufficiently 
smooth over the receiver bandwidth, then we can write 

NL = NSL +10 log w . 

Otherwise, the total bandwidth must be broken down into smaller 
bands, the band level found for each, and then these band levels 
combined according to the nomogram for combining independent 
levels. 
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REVERBERATION 



If the sound field generated by the source is strong enough, and 
there are enough scatterers present in the vicinity of the desired 
target, then unwanted reflections and scattering of sound from 
the undesired scatterers are received along with the echo from 
the desired target. The echo level of the signal received from 
the desired target is 

EL * SL - 2TL + TS (monostatic) 

where TS is the target strength of the desired target. The 
reverberation level RL can be written in the form 

RL * SL - 2TL + TS’ (monostatic) 



where TS’ is the combined target strength of all the undesired 
targets which are generating the masking level. Since the RL 
is the masking level when reverberation is the dominant source 
of unwanted sound, the SONAR equation becomes 



or 



EL - RL ^ DT 
TS - TS’ ^ DT 



(monostatic ) 



The central problem in obtaining TS’ is determining the volume 
or surface area containing the scatterers which is illuminated 
by the source, sending sound to the receiver, and sending the 
sound so that it arrives at the same time as the signal from the 
desired target. This is suggested in the sketch below. 
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VOLUME REVERBERATION 




The target strength TS* of the undesired scatterers can be 
written as 

TS’ * 10 log s y V 

where V is the volume within which the reverberation sent to 
the receiver is in competition with the reflected signal from 
the desired target and s v is a measure of the scattering cross 
section per unit volume. It has become conventional to express 
s v in dB terms through the definition 

S v = 10 log s v 

where S v is the scattering strength . Thus, we have 

TS’ * S v + 10 log V. 

The calculation of the reverberating volume V contains the heart 
of the problem. As suggested previously, this is a study in 
geometry, involving the way in which major lobes of source and 
receiver overlap. 

If we restrict ourselves to the monostatic situation, and assume 
that the major lobe of the receiver completely overlaps the 
major lobe of the source, it then follows that the volume V is 
given by 

10 log V * 10 log + TLg + to log fl 

where t is the time duration of the acoustical pulse and lb is 
the solid angle subtended by the major lobe of the radiation 
pattern of the source. Recall that this is related to the 
directivity index by DI ~ 10 log(47T/n.). The term TLg arises 
from the fact that the cross sectional area of the reverberating 
volume clearly depends on the range from the source, and this can 
be seen to be expressed by the geometrical spreading of the sound 
beam. Thus, for this simple monostatic situation, we have 

TS’ - S v + 10 log |t + TLg + io log 4 tt - DI 
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SURFACE REVERBERATION 




In the case of scattering which arises from a surface (such as 
the sea surface or the bottom) or from some well-defined layer 
such as the bubble layer underlying the water-air interface, 
the reverberation target strength is expressed as 

TS* = S g + 10 log A 

where A is the surface area generating reverberation in competition 
with the desired target and S s is the scattering strength (per 
unit area). 

As before, the problem is one primarily of geometry, involving 
the way in which source and receiver major lobes overlap. If we 
again restrict ourselves to the simple monostatic geometry used 
previously (for which the major lobe of the receiver completely 
overlaps that of the source), then the surface area of interest 
can be seen to be 



A = i r 9 c t 

where 9 is the angle in radians of the horizontal angular width 
of the radiation pattern of the source. This then gives the 
equation 

TS* = S o +10 log £t + io log r + 10 log 9 
s 2 
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REPRESENTATIVE SCATTERING STRENGTHS 



The values quoted here are representative only. Large 
deviations can occur depending on location, season, biological 
activity, and so forth. 

Deep Scattering Layers? The scattering strength S v 
can vary between about -90 to -60 dB within the 
various sublayers. 

Water Volumes In the clear open ocean waters, S y 
can be found to range between -100 to -?0 dB. 

Surface Scatterings For reasonably low grazing angles, 
the combined scattering from the ocean surface and 
the population of bubbles found at shallow depths 
below the surface gives scattering strengths S s 
ranging around -50 to -30 dB for sea states 
between 1 and 4. 

Bottom Scatterings For grazing angles between 20° 

and 60° , representative values of the scattering 
strength S s range between -40 and -10 dB. The 
exact values depend very strongly on the particular 
bottom and its composition. In general, however, 
as the grazing angle goes to zero, the scattering 
becomes negligible and S g goes to - oo dB. 



REVERBERATION LEVEL DEPENDENCIES 

Notice that for either volume or surface reverberation 
the reverberation target strength TS’ increases with increasing 
range, either as TL„ (volume) or 10 log r (surface). This is to 
be compared with ambient noise masking, for which there is no 
range dependence. This means that, given sufficient source level, 
as -SL 'MincreasesSONAR performance will change from ambient - 
noise limited performance to reverberation limited performance . 
While detection for ambient noise limited performance depends 
on the source level and directivity index, notice that when the 
performance becomes reverberation limited there is no dependence 
on the source levels additional power will not improve 
reverberation limited performance. Reverberation can be reduced 
by increasing the directivity of the source and/or receiver, or 
by decreasing the pulse length, given by ct. Increased directivity 
has the disadvantage of reducing scanning rates, and decreased 
pulse length may result in reduction of the target strength TS 
of the desired target. 



104 



DOPPLER SHIFT FOR ACTIVE SONAR 



In the case of active SONAR, any motion of the platform 
containing the source will cause the signal propagated into 
the water to have a frequency shift which depends on the 
bearing into which the signal is sent. This doppler-shifted 
signal is subjected to a further doppler-shift when it is 
reflected from the moving target. Without going into the specifics 
of the various frequency shifts, we shall simply quote some results 
of interest: 

Because of the two-way propagation of acoustical energy from 
the moving source to the moving target and then reflection back 
to the moving receiver, the required bandwidth for the receiver 
of an active sonar operated under monostatic conditions is 

w in Hz = 1.4* (Range rate in knots )• (Frequency in kHz) 

Notice that this is twice the bandwidth required in the passive 
SONAR case. 

If the target is stationary with respect to the surrounding 
water, then the desired echo will have the same frequency as the 
reverberation and the ability of the receiving system to recover 
signal from noise will be minimized. On the other hand, when the 
target is moving through the water, then the echo and reverberation 
will have different frequencies and detection will be much easier. 
This effect is quite easy to observe by ear: The detectibility 

of an "up-doppler" or "down-doppler" echo in reverberation is 
very easy even when the relative loudness of the desired echo is 
rather weak. When there is no frequency change between echo and 
reverberation, however, the ability of the ear to detect the 
echo is drastically reduced. The strong increase of the detection 
threshold for slow targets points to the use of parallel filters 
of relatively narrow bandwidths to facilitate the separation of 
signal from reverberation when the difference in frequency is 
small. 
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SL 



2TL + TS ML + 



DETECTION THRESHOLD 



Energy Detection 

If, as was the case for the passive SONAR equations, the 
mechanism for the detection of the desired signal is based on 
detecting an excess of acoustical energy in the frequency band 
of interest, then we have the same expression for the detection 
threshold , 



DT = 5 log 

where t is the integration time. In the case of active sonar, 
t should be the length of the received echo. If the integration 
time is less than or greater than this pulse duration, it can be 
shown that the DT is increased above this expression, and 
performance degraded. 

Correlation Detection 



An alternative mode of signal processing is possible in 
the case of active SONAR. Since the amplitude and frequency 
properties of the tone burst generated by the source are known, 
it is possible to search for a signal of these same properties 
in the received echo. If. the detailed shape of the received 
echo matches that of the sent pulse, then it can be shown that 
the detection threshold is given by 

DT = 10 log ^ 

when the processing is that of correlation of the received signal 
against a replica of the generated pulse. There are many 
modifications of this basic idea, but all rely on the technique 
of multiplying the received signal by a time-delayed model of 
the sent pulse and integrating the resultant product over the 
pulse duration. This is done for monotonically increasing 
delay times, and the resulting function of delay time studied 
for the presence of a sharp peak representing the matching of 
the delayed replica with the same signal appearing in the received 
signal and noise. In practice, the finite size of the target, 
multipath interference, and the fluctuations introduced by the 
inhomogeneities of the transmission of the signal through the 
water alter the received echo so that the correlation process 
is not as good as would be desired. As a result, the DT for 
a real system operating in the real world is more than the 
theoretical expression would predict. 
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SUMMARY OF THE ACTIVE SONAR EQUATIONS 



Ambient noise limited conditions 



SL - 2TL + TS ^ H5L - AG + 10 log w + DT 



Reverberation limited conditions 



ts - ts' ^ nr 
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EXAMPLE 



An active SONAR operates at 1 kHz with a source level of 220 dB re 1 fJ Pa, 
a directivity index of 20 dB, a horizontal beam width of 10°, and a pulse 
length of 0.1 sec. Correlation detection is used and it is desired to have 
P(D) - 0.50 with P(FA) * 10"^. The target strength of the submarine is 
expected to be 30 dB and its speed may be up to 20 kts. Both the SONAR 
and the submarine are in a mixed layer of depth 100 m with the source at 
36 m. The sea state is 3> the scattering strength for surface scattrrers 
is -30 dB, and volume scattering is negligible. For simplicity, assume that 
conditions are such that the detection threshold is the same for both noise- 
limited and reverberation-limited performance. 

Solution: 

w ■ l.U (frequency in kHz) (range rate in kts) ■ 1.1* x 1 x 20 - 28 Hz 

d » 15 same calculation as in example 1 (passive) 

DT » 10 log ( d/2wt ) - 10 log (15/2x28x0.1) - 10 log 2.68 - 5 dB 

NL - NSL - DI + 10 log w - 62 -20 + Uj. - 56 dB 

r t ■ 1310 m and r g ■ 84 OO m same calculation as Example 3 (passive 

Assume that the detection is noise limited 

2TL ^ SL + TS - NL - DT - 220 + 30 -56 _ 5 - 189 

TL ^ 95 

-4 

TL ■ 31 + 10 log r +2.6 x 10 r same calculation as Example 3 

10 log r + 2.6 x 10 ^ r = 95-31 = 64 



r 


10 log r 


2.6 x 10” 4 r 


TL 


10 4 


40 


26 


43 


10 5 


50 


26 


76 


5 x 10 4 


47 


13 


60 


6 x 10 4 


48 


16 


64 



r = 6 x 10 4 m if noise limited 
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EXAMPLE (CONTINUED) 



But, what is the RL at this range? 

TS' - S g + 10 log ( ct/2 ) + 10 log r + 10 log 6 

- -30 + 10 log (1500 x 0.1/2) + 48 + 10 log 0.171* 

“ -30 +19+48-8 * 29 dB 

RL - SL - 2TL + TS' 

» 220 - 189 + 29 . 60 dfl 

So 

RL ^ NL 

Therefore the detection will be reverberation limited 
TS - TS » DT 

TS » - -30 + 19 + 10 log r - 8 - -19 + 10 log r 

TS - TS' 'V DT 
30+19 - 10 log r >, 5 
10 log r ^ 14* 

r - 2.5 x lO* 1 m 
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APPENDIX A: LOGS 



All you have ever wanted to know about the decibel but were afraid to ask . 
The language in which the SONAR equation is written is the language of the 
} decibel. This self-paced exercise is intended to make you familiar with the 

grammar of this language. 

Instructions . For each numbered section, first attempt to do the problems. 

(The answers are on the last page.) If you have no troubles, skip over to the 
next section. On the other hand, if you have difficulties^ read through the 
discussion and then try the problems again. Go to the next section only when 
you can do the problems without difficulty. By the time you finish all the 
sections, you should be speaking "decibel" like a native. 

1. EXPONENTIAL NOTATION 

Problems A 

1. Write each of the following numbers ir. exponential forms 

a) 2,365,000 

b) 872 

c) 3340 

d) 42.1 

e) 600 

f) 1.0 

2. Write each of the following numbers in conventional form 

a) 2.689 x 10 2 

b) 7.30 x 10 2 

c) 6.29 x 10 4 

d) 8.9 x 10 6 

e) 3.261 x 10 2 

f) 1.00 x 10° 
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3. Write each of the following numbers in exponential form 

a) 0.00252 

b) 0.012 

c) 0.00002 

d) 0.01002 

e) 0.00328 

4. Write each of the following numbers in conventional form 

a) 1.03 x 10" 6 

b) 3.2 x 10~ 2 

c) 6.89 x 10 -4 

d) 1.003 x 10 _1 

e) 1.00 x 10 _1 



A succinct way of expressing very large and very small numbers is by 
exponential (or scientific) notation. For example, one-hundred and twenty-six 
million can be written either as 

126,000,000 

or as 

1.26 x 10 8 , 

where the exponent 8 denotes the number of places that the decimal must be 
moved to the right to express the number in its conventional form. 

There is no unique way of expressing a number in exponential form, since 
1.26 x 10® = 12.6 x 10 7 = 0.126 x 10 9 = etc. 

Which form you use depends on convenience. 

Numbers smaller than unity can be handled in a similar manner. For 
example, one-hundred and twenty millionths can be written either as 

0.000120 

or as 

' 1.2 x 10‘ 4 , 
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where the negative exponent denotes that the decimal place is to be moved to 



the left . 

A special case that must be watched is 

1 = 1 x 10° = 10° 

which is consistent with the above conventions. 
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2 . ALGEBRA OF EXPONENTS 



Problems B 

1. Calculate the following and express the answer in exponential form 

a) 6.15 x 10 3 + 2.34 x 10 3 

b) 3.2 x 10 2 - 1.46 x 10 3 

c) 9.1 x 10 4 + 2.4 x 10 5 

d) 3.58 x 10" 2 + 1.26 x 10" 2 

e) 6 x 10" 6 + 3.281 x l(f 4 

f) 3.5 x 10 4 - 2 x 10 1 

2. Calculate the following and express the answers in exponential form 



a) 


2.50 


X 


10 6 x 2.00 x 10 4 


b) 


4.40 


X 


10 4 / 2.20 x lo 5 


c) 


3.26 


X 


10 3 x 2.00 x lO -2 


d) 


5.34 


X 


10 -2 / 2.67 x 10 5 


e) 


8.68 


X 


10 _1 / 4.34 x 10~ 2 



3. Calculate the following and express the answer in exponential form 

a) (1.1 x 10 4 ) 2 

b) (2.0 x 10” 3 ) 4 

c) (1.44 x 10 2 ) 1/2 

-2 1/2 

d) (1.44 x 10 ) ' 

e) (8.0 x 10 6 ) 1/3 

f) (2.56 x 10 4 ) 0 * 5 

. . -6. 0.3333 

g) (8 x 10 ) 

To add and subtract numbers expressed in exponential form it is necessary 

4 

to express both numbers in the same exponent. For example, adding 3.56 x 10 

3 

to 2.1 x 10 is accomplished by 

3.56 x 10 4 
0.21 x 10 4 

3.77 x 10 4 
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When asked to subtract 8.67 x 10 1 from 2.23 x 10 3 , note that the answer is 
2.23 x 10 3 , for the term to be subtracted affects the larger term beyond the 
figures retained. 

Exponential notation really shows its worth in calculations involving 
multiplication and division. We see that 

120 x 2000 = 240,000 

or 1.2 X 10 2 X 2 x 10 3 = 2.4 x 10 5 

and 

36000 / 200 = 180 
or 

3.64 x 10 4 / 2 x 10 2 = 1.8 x 10 2 . 

The rules are simple: 

In multiplication, exponents add. 

In division, exponents subtract. 

A special case that is often useful is 

2 0 2 -2 
1 / 10 = 10 / 10 ^ = 10 , 

which allows exponents to be moved between numerator and denominator. 

Power and roots are trivial in exponential notation: 

(1.2 x 10 3 ) 2 = 1.2 x 10 3 x 1.2 x 10 3 = 1.44 x 10 6 , 



or for any power n 



and. 



or for any power of n 



doV = 10™ 



(1.44 X 10 6 ) 1/2 = <1.44> 1/2 X (10 6 ) 1/2 ' 2 ' 2 X 10 ' 



do 1 ") 1 /" - io n/n . 



This latter form introduces the possibility of fractional exponents; these will 
be discussed in Section 4. 



1 

l 
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3. LOGS (IN BASE TEN) 

Problems C 



Find x 




a) 


log 


10 3 = x 
-6 


b) 


log 


10 = x 


c) 


log 


10 X = 4 


d) 


log 


10 X = -3 


e) 


log 


io x = 0 



Since we say that "ten to the fourth power equals ten thousand" we can 
ask "to what power must ten be raised to get ten thousand?" The answer is 
obviously "four." You are now an expert on logarithms, since the "log" of 
10,000 is 4. This is written 

log 10,000 = 4 or log 10 4 = 4 
and this is all there is to logs . 

Find the log of x by asking yourself to what power must I raise 10 to 
get x. 

It is equally appropriate to ask "the log of what is 4," that is what 
x 

is x so that log 10 =4. 
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4. DECIMAL EXPONENTS 



PROBLEMS D 

1. Find 

a) log 4 

b) log 40 

c) log 0.4 

d) log 4 2 

e) log 1/4 

2. Find the number whose log is 

a) 0.70 

b) 1.70 

c) -1.30 

d) -1.70 

e) -2.70 

So far we have assumed integer exponents, but everything we have 
said also is true for decimal exponents. For example, 

10 °' 30 = 2.0 

and 



log 2 - 0.30 

The table below will allow you to calculate, with a fair degree of accuracy, 
all logs. 



X 


1.0 


1.2 


1.5 


1.7 


2.0 


.3.0 


4.0 


5.0 


6.0 


7.0 


8.0 


9.0 


10.0 


lo<£ X 


0.00 


0.08 


0.18 


0.23 


0.30 


0.48 


0.60 


0.70 


0.78 


0.85 


0.90 


0.95 


l.Oi 
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If you ever need more accuracy you will have to use slide rule, calculator, 
or tables. 

To find the log of a number greater than 10 proceed as follows: 
log 7000 = log 7 x 10 3 = log 7 + log 10 3 
= 0.85 +3 = 3.85 

or 

log 0.003= log 3 x 10 “*= log 3 + log 10 ^ 

= 0.48 - 3 = -2.52 

There is one property of logs that, if used correctly, will greatly 

simplify many calculations. This will be illustrated by an example: 

2 

log 36 = log 6 = log (6x6) = log 6 + log 6 = 2 log 6= 1.56 

( Calculate log 36 and compare to 2 log 6.) 

In general 

i n 

log x = n xog x. 

It is also true that 

log xy = log x + log y. 

log 36 = log (4x9) = log 4+2 log 3 = 0.60 + 0.96 = 1.56 
But , remember that 

log (x + y) is not equal to log x + log y ! 1 I 



117 



5. THE DECIBEL 



PROBLEMS E 



1. Express these signal- to-noise ratios in dB 

a) 2 

b) 1/2 

c) 10 

d) 1/10 

e) 100 

2. The voltage of a signal is 1 volt*, express the signal level with 
reference to the following voltages 

a) 10 -12 

b) 10" 6 

c) 1 

d) 10 6 



By convention, when the log of a pressure or voltage is multiplied 
by 20, the resulting number is said to be measured in "decibels” , donated 
by dB. 

20 log 2 = 6dB 
20 log 10 = 20dB 

Decibels are used in two different ways: 

a) to express a ratio such as the signal-to-noise ratio 

signal- to-noise ratio = 20 log (signal voltage/noise voltage) 

b) to express a level such as signal level 

signal level = 20 log (signal voltage/reference voltage) 

In the latter case it is absolutely essential to specify the reference, 
e.g., "the signal level is 120dB re 10 ^ volts". 
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Problems A: 



Problems B : 



Problems C: 



Problems D: 



Problems E: 



ANSWERS 

1. a) 2.365 x 10 6 , b) 8.72 x 10 2 , c) 3.34 x 10 3 , d) 4.21 x 10 1 , 

e) 6 x 10 2 , f) 1.0 x 10° 

2. a) 268.9, b) 730, c) 62,900, d) 8,900,000, e) 326.1 

f) 1.00 

3. a) 2.52 x 10 -3 , b) 1.2 x 10 -2 , c) 2 x 10~ 5 , d) 1.002 x 10~ 2 
e) 3.28 x 10 -3 . 

4. a) 0.00000103, b) 0.032, c) 0.000689, d) 0.1003, e) 0.100 



l.a) 


8.49 x 10 3 , b) -1.14 x 10 3 , 


c) 3.3 x 10 5 , d) 4.84 


x 10 -2 


e) 


3.34 x 10 -4 , f) 3.5 x 10 4 . 






2. a) 


5.00 x 10 10 , b) 2.00 10 _1 , 


c) 6.52 x 10 1 , d) 2.00 


x 10 -7 , 


e) 


2.00 x 10. 






3. a) 


1.2 x 10 8 , b) 1.6 x 10 -11 , 


c) 1.2 x 10 1 , d) 1.2 x 


io" 1 . 


e) 


2.0 x 10 2 , f) 1.6 x 10 2 , g) 


2. x 10 -2 . 





l.a) 3, b) -6 ,c) 4, d) -3, e) 0 

1. a) 0.60, b) 1.60, c) -0.40, d) 1.20, e) -0.60. 

2. a) 5.0, b) 50, c) 0.050, d) 0.020, e)0.0020 

1. a) 6dB , b) -6dB, c) 20dB, d) -20dB, e) 40dB. 

2. a) 240dB , b) 120dB, c) OdB, d) -120dB, e) -240dB. 
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